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We evaluate the ratio of the shear viscosity to the relaxation time of the shear flux above but 
near the critical temperature T c in SU(3) gauge theory on the lattice. The ratio is related to Kubo's 
canonical correlation of the energy-momentum tensor in Euclidean space with the relaxation time 
approximation and an appropriate regularization. Using this relation, the ratio is evaluated by direct 
measurements of the Euclidean observables on the lattice. We obtained the ratio with reasonable 
statistics for the range of temperature 1.3T C < T < 4T C . We also found that the characteristic speed 
of the transverse plane wave in gluon media is almost constant, v ~ 0.5, for T > 1.5T C , which is 
compatible with the causality in the second order dissipative hydrodynamics. 
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I. INTRODUCTION 

One of the highlights at RHIC (the Relativistic Heavy 
Ion Collider) is a success of relativistic hydrodynamic 
models for the description of the space-time evolution of 
the quark-gluon matter created in heavy ion collisions 
P-Q. This novel feature revealed that the hot quark- 
gluon matter is a strongly correlated system, which is 
one of the most significant subjects over broad areas in 
physics. Heavy ion collisions at LHC (the Large Hadron 
Collider) have already reported new experimental results. 
The transverse momentum dependence of the elliptic flow 
observed at LHC is similar to that at RHIC [4[ . The ex- 
perimental result implies that the hydrodynamic descrip- 
tion for the quark-gluon matter is as good at LHC as at 
RHIC. 

Recently the importance of dissipative effects in the 
hydrodynamic models has been recognized 5, 6]. The 
simplest relativistic dissipative hydrodynamics is the first 
order theory, which in the non-relativistic limit reduces 
to Navier- Stokes theory @, The first order theory is, 
however, accompanied with the acausal problem and in- 
stability in numerical simulations [9|. One of the strate- 
gies to evade these problems is to extend the theory to 
the second order with regard to dissipative fluxes. In 
the second order theory, however, there appear many 
phenomenological parameters as second order transport 
coefficients in addition to the first order ones. These 
transport coefficients cannot be determined within the 
framework of hydrodynamics. Ab initio calculation based 
on microscopic theory, i.e. Quantum Chromodynamics 
(QCD) in our case, is required to study and constrain 
the parameters in dissipative hydrodynamic models. 

Since the temperature range realized at RHIC and 
LHC is not within the reach of perturbative QCD, we 
need to employ non-perturbative approaches to investi- 
gate the transport coefficients. At present, lattice sim- 
ulation is the only systematic way to calculate physical 
quantities in such a non-perturbative region. There are 
several pioneering studies which analyzed transport co- 
efficients on the lattice [T5Mx3 | . These studies have used 
Green-Kubo formulas, which relate the transport coeffi- 



cients to the low energy behavior of the corresponding 
spectral functions. In this method, one needs to ex- 
tract the spectral functions from Euclidean correlators 
obtained on the lattice. This step is, however, nontrivial 
because it is an ill-posed problem [15j . 

In this paper, we focus on the ratio of the shear 
viscosity to the relaxation time of the shear flux. In 
Refs. [Ill [13] , it is argued that the ratio is related to 
a static fluctuation of the energy-momentum tensor by 
rewriting Green-Kubo formula for a classical system with 
the relaxation time approximation. We modify their ar- 
guments to treat a quantum system and show that the 
ratio is related to Kubo's canonical correlation of the 
energy-momentum tensor. It enables us to evade the dif- 
ficulty in the analysis of spectral functions and calculate 
the ratio on the lattice directly. Whereas the ratio itself 
is not a transport coefficient, its determination reduces 
the number of the free phenomenological parameters in 
hydrodynamic models. We note that a similar argument 
has been also given in Ref. [l8j . 

We need to regularize Kubo's canonical correlation to 
extract physical quantities, because it is ultraviolet di- 
vergent. We perform the regularization by taking the 
vacuum subtraction. We argue that one further needs to 
take care of the contribution arising from contact terms 
to obtain physically meaningful quantities in addition to 
this prescription. While Green-Kubo formula does not 
contain the equal space-time operator by definition, the 
lattice observable includes its contribution. This con- 
tribution is usually removed solely by the vacuum sub- 
traction, because the coefficients of the contact terms 
are not dependent on the temperature. For the case of 
the energy-momentum tensor, to the contrary, the coef- 
ficients of the contact terms in Kubo's canonical correla- 
tion of the energy-momentum tensor are proportional to 
the energy-momentum tensor themselves, and their ex- 
pectation values are dependent on the temperature. We 
thus have to remove the contribution of the contact terms 
separately besides the vacuum subtraction. 

The structure of this paper is as follows. In Sec. [TT1 we 
review Israel-Stewart theory [l9j] , which is one of the sec- 
ond order dissipative hydrodynamics. We show how to 
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introduce the relaxation times of the dissipative fluxes to 
hydrodynamics. In Sec. lIIII we relate the ratio of the shear 
viscosity to the relaxation time to Kubo's canonical cor- 
relation of a spatial component of the energy-momentum 
tensor. We discuss how to remove the unphysical diver- 
gence in the correlator and to extract the ratio from the 
canonical correlation in Sec. IIV1 In Sec. |V]we formulate 
a method to calculate the ratio on the lattice and show 
numerical results in SU(3) lattice gauge theory. The last 
section is devoted to conclusions and discussions. 



II. ISRAEL-STEWART THEORY AND 
RELAXATION TIME 

In this section, we give a brief review on Israel-Stewart 
(IS) theory and the role of the relaxation times in this 
theory. Hydrodynamics consists of equations of motion 
for conserved currents. The only conserved current in 
pure gauge theory is the energy-momentum tensor T^, 
for which the conservation law reads 



= 0. 



The energy-momentum tensor is decomposed as 
= eu^u" - (P + II)( 5 ^ - u^u") + TTfn 



(1) 



(2) 



where e and P are the energy density and pressure in 
equilibrium, and II and are the dissipative fluxes for 
the bulk and shear channels, respectively. The metric 
g» v is defined by g» v = diag(+, -, -, -). In Eq. ©, 
we took the energy frame [8|, i.e. is proportional to 
energy-flow with the normalization u^u^ = 1. 

In second order theory, II and are regarded as 
dynamical variables. One of the arguments to obtain 
evolution equations for these variables besides Eq. (fTJ) is 
a phenomenological one [l9l - [2"lj |. In the derivation, the 
evolution equations are obtained as a consequence of the 
second law of thermodynamics 



>0, 



(3) 



where — sw M and s is the entropy density in the energy 
frame. By choosing the entropy density to be the equilib- 



rium one, s 



'cqj 



we obtain the first order theory with 
Eq. ([3]) and linearity between the dissipative fluxes and 
thermodynamic forces. To extend the theory to second 
order, it is assumed that the entropy near equilibrium is 
given by a function of II and Ttuu, and that it is analytic 
in equilibrium II = tt^ — [19j . The entropy density up 
to the second order in Taylor expansion with respect to 
II and 71^,1, then reads 



2T 



(/3 n 2 + fan^) 



(4) 



where /3o and 02 are phenomenological parameters. The 
second law of thermodynamics Eq. ([3]) with Eq. (0| leads 
to a set of evolution equations, which are referred to as 



IS equations. In the shear channel, for example, the evo- 
lution equation reads 



1 r i 
TTuk = Wuu + Oi\ , 



(5) 



where ■h liv = dn^fy/dt and a denotes contributions 
which vanishes in a uniform medium. Eq. ([5]) shows that 
7r MI/ damps with a time scale 7v in a uniform medium, t v 
therefore is called the relaxation time of the shear flux 
and is related to f3 2 in Eq. Q as 



(6) 



where r\ is the shear viscosity }19| . The corresponding 
first order equation is obtained by taking the limit r n — > 
in Eq. Similarly the evolution equation in the bulk 
channel includes the relaxation time of the bulk flux tjj, 
which is related to 0o in Eq. ((4]) as rn = CA)- The second 
order theory can become causal due to these relaxation 
times of dissipative fluxes. 

The evolution equation Eq. ([5]) implies that the tem- 
poral correlator of the shear flux behaves as 



(7fp„(*)7f^(0)) ~ e- t/T " (7^(0)7^(0)), 



(7) 



where tr^i) = y f v d 3 x Tr^(t, x) is the spatial average 
of 7r^„ in a volume V. We note that TT^ v {t) in Eq. (J7J is a 
classical variable in hydrodynamical equations. We will 
later argue an extension of this relation to the quantum 
correlator. 

In addition, we introduce the wave front speed (char- 
acteristic speed) for the transverse plane wave v defined 
by M 



V 



V 



r w (e + P) tvTs 
For the causality to be satisfied, v 2 < 1 is required. 



(8) 



III. RATIO OF THE SHEAR VISCOSITY TO 
THE RELAXATION TIME 

A. The ratio r?/r w 

We now relate the ratio rj/r^ to Kubo's canonical cor- 
relation of the energy-momentum tensor in Euclidean 
space. The original idea of the following discussion was 
proposed within the framework of classical theory in 
Refs. [lfl [I?} • We modify their arguments to treat quan- 
tum systems. 

We start from Green-Kubo formula in the shear chan- 
nel 



T 1 =-Jd 3 xJ dtt([TT 12 (t,x),TT 12 {0,0)}) 



-V 



dtt <[7fl 2 (t),7fl 2 (0)]>, 



(9) 
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where (O) = Tr[e _ ' ff / T 0]/Tr[e _if/T ] is the statistical av- 
erage in equilibrium with H being the hamiltonian. Here 
we have set Boltzmann constant fee = 1. To rewrite this 
formula, we use the following identities, 



[e- H/T ,A] = 
([A,B])=Tr 



■" 1 l 1/T d\e XH \A,H}e- XH , 



(10) 



,-H/T 



[A,B] 



h/T 



dXA(-iX)B). 

(11) 

Substituting Eqs. (TO) and dH} to Eq. ©, we obtain 



V 



h/T 



V = ~jJ dtt(j dA7ria(t-iA)7fi 2 (0)) 



o 



d\ 7ria(*-tA)7fi3(0)). 



(12) 



Eq. (IT21) relates the shear viscosity to an integral of 
the correlator over the complex time t — iX. In order to 
modify this relation further in a way feasible in lattice 
calculations, we assume that the imaginary time integral 
of the correlator damps exponentially with real time, 



dX 7fi2(t-iA)7fi 2 (0)> 
e-*/ T "( / dA7fi 2 HA)7fi 2 (0)) r 



(13) 



for t > 0. On the right hand side of Eq. (|T3|) . we have put 
the subscript "reg" to the statistical average, because the 
expectation value (J dA7fi 2 (— iA)-7f 12 (0)) is in fact ultravi- 
olet divergent and ill-defined, and should be regularized. 
This expectation value should not contain the contribu- 
tion arising from the contact term, because the original 
formula Eq. ^ does not include it. We postpone detailed 
discussion on these issues to the next section, and for the 
moment assume that the right hand side of Eq. (|13l) can 
be defined appropriately. 

In the following, we refer to Eq. (fT3|) as the relaxation 
time approximation, owing to its analogy to the standard 
approximation in nonequilibrium statistical mechanics, 
(O(i)O(O)) ~ e-*/ r (Cl(0) 2 ). Note, however, that Eq. (13} 
is slightly different from this approximation; because in 
Eq. ([TB"} the correlator between operators separated along 
complex time is concerned, while the standard relaxation 
time approximation is for a real time separation. Since 
Eq. (TIB")) should reduce to Eq. ([7]) for the range of t where 
the second order hydrodynamics is applicable, the time 
scale appearing in the exponential function in Eq. (fT3"|) is 
identified to be the relaxation time of the shear flux tv 
in Eq. (5}. 

Using Eq. (fl3"}. Eq. ([12} reduces to 

V r°° r h / T 
V-jrJ dtj dA e -*/ r ^ 12 HA)^ 12 (0)) rcg 

v r H/T 

rfA(7f 12 HA)7fi 2 (0)) reg . (14) 



In the local rest frame of a medium, one can replace 
7r 12 with a spatial component of the energy-momentum 
tensor Ti 2 . Eq. (|14[) then means that the ratio T)/t v is 
given by 



H 

T-rr 



V 
H 



h/T 



dX <f 12 (-*A)T ia (0)) 



reg 



T 



(^i 2 )reg- 



Here (A 2 ) denotes Kubo's canonical correlation 



h/T 



dX(A(X)A(0)}. 



(15) 



(16) 



The same argument applies to the bulk channel, which 
leads to a formula to represent the ratio C/ T n with di- 
agonal components of the energy-momentum tensor. We 
note that Eq. (IT4|) is also derived by the projection oper- 
ator method 1181. 



B. Classical limit 



In the classical limit H/T 

integral with respect to A in Eq 
Eq. (fT5|) then reduces to 



0, one can replace the 
J Q h/T dX by h/T. 



= ^(T 1 2 2 ) 



(17) 



7;/tv is, therefore, given by a static fluctuation of Ti 2 in 
classical theory. Eq. (fTTl) i s the formula derived in the 
original works Refs. [Iq.ll7|. 

In a classical system, one can also derive Eq. (fTTj) us- 
ing Einstein principle and the entropy density in the IS 
theory. According to Einstein principle, the probability 
distribution P(x) of a state variable x in equilibrium is 
given by P{x) ~ e styX ^ where S(x) is the nonequilibrium 
entropy in a volume V. If one identifies the IS entropy 
Eq. (01 to be that in Einstein principle, the probability 
distribution for a state variable 7fi 2 in equilibrium state 
in a volume V is given by 



PAnw) ~ exp 



(18) 



With this distribution function the fluctuation of 7r 12 is 
calculated to be 



'12/ 



/ dTfnj 2 P^(Tr) T rjT 
J dfiPjf) ~ 2f3 2 V 



(19) 



which is equivalent with Eq. (I17p . It should, however, be 
remembered that Einstein principle is applicable only to 
classical systems. In quantum mechanics, when a state 
variable x does not commute with the hamiltonian, x and 
the energy cannot be determined simultaneously. This 
means that the entropy cannot be defined as a func- 
tion of non-conserving state variables in quantum sys- 
tems (22)]. Einstein principle therefore is not applicable 
for such cases. 
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IV. REGULARIZATION 

From now on we set K = 1. On the lattice, rj/r^ 
is obtained by direct measurement of Kubo's canonical 
correlation of the spatial off-diagonal component of the 
energy-momentum tensor 

(T^) = yJd i x{T 12 (x)T 12 (0)}, (20) 

where x — (t,x). This canonical correlation is, however, 
different from (T 1 2 2 ) rog in Eq. ([T5j) because of two rea- 
sons. First, the canonical correlation Eq. (|2"0|) diverges 
in the ultraviolet limit. This divergence originates from 
the short distance behavior of the Euclidean correlator of 
Ti2(x). Second, Eq. (|20|) contains the equal space-time 
correlator, i.e. the contact term, while Eq. (|15p does not. 
The absence of the contact term in Eq. (fl~5j) is a conse- 
quence of the fact that Green-Kubo formula Eq. ((9|) does 
not include the equal space-time due to the commutator. 
To evaluate tj/t^ on the lattice with the measurements 
of Eq. (l20l) . one has to take account of these differences. 

Since both of these effects originate from the short dis- 
tance behavior of the Euclidean correlator of Ti2(x), we 
examine them by the operator product expansion (OPE). 
According to the OPE calculation [ijjj], the correlator at 
short distance behaves as 

(T 12 (x)T 12 (0)) * C T ^ + ]T C,(T^(0))5^(x) + • • ■ , 



where Ct and C M are c-number Wilson coefficients de- 
termined by perturbation theory. Because of the dimen- 
sional reason, the first term on the right hand side in 
Eq. (f2~Tj) is proportional to |x|~ 8 . This term leads to 
the ultraviolet divergence of Eq. (j2T)]) . This term neither 
has medium effects nor affects long distance behavior re- 
sponsible for hydrodynamics. Therefore, this divergence 
should be eliminated when one discusses transport coef- 
ficients. The contribution is removed by subtracting the 
correlator in the vacuum, 

J d 4 x (T 12 (x)T 12 (0)) 

= J d 4 x ((T 12 (x)T 12 (0)) T - (T la (ar)T 12 (0))r=o)- (22) 

The next terms on the right hand side of Eq. (f2~Tj) are 
proportional to dimension four operators. These terms 
are the contact terms, i.e. they are proportional to 
S^'(x), as presented in Eq. ([2~Tj) . As discussed already, 
we need to remove these terms from Eq. (1201) when eval- 
uating Eq. (fTS"]) . i.e. 

j d 4 x (T 12 (x)T 12 (0)) rcg 

= J d 4 x ({T 12 (x)T 12 {0)) - i(contact terms)^) . (23) 



TABLE I: Simulation parameters. N a and N T are the num- 
bers of the lattice sites in the spatial and temporal directions, 
respectively, a and L a are the lattice spacing and spatial 
lattice size, respectively. 



P = 6/<7 2 


N a 


N T 


a [fin] 


La [fm] 


6.499 


32 


4,6,8,32 


0.046 


1.5 


6.205 


32 


4,6,8,32 


0.068 


2.2 


6.000 


32 


4,6,8,16 


0.094 


3.0 


6.000 


16 


4,6,8, 16 


0.094 


1.5 



The statistical average of the contact terms has tempera- 
ture dependence, because it is proportional to the energy- 
momentum tensor. This contribution therefore cannot be 
removed by the vacuum subtraction Eq. (|2"2"j) . In Eq. (|23l) . 
we subtracted the contribution of the contact terms with 
a factor 1 /2 because of the following reason. As shown in 
Eq. (|2T|) . the support of the contact terms or (^-functions 
is located at the origin. The integral interval for the 
time in Green-Kubo formula Eq. © is < t < oo, and 
the same holds for Eq. ([20]) as well. Thus f]/r n contains 
the half-contribution from the vicinity of the origin. The 
higher order terms in Eq. (I21[) shown by the dots neither 
yield divergence nor are singular. 

In the leading order OPE calculation 0, [25| , the con- 
tact terms are most conveniently evaluated by taking the 
high frequency limit of the correlator with zero spatial 
momentum 

G (q 4 ->■ oo, q = 0)=J2 C^(T^(0)) Q 

A" 

= ^oo)o + ^ 2 )o, (24) 

where G Q (q A ,q) = J d 4 x e l ^ T -^ s) (T 12 {x)T 12 (0)) and 
(F 2 ) = {F^ V F^ V ) is the gluon condensate. In the first 
equality of Eq. (|24[) , the higher order terms in Eq. (f2Tj) 
vanish by Riemann-Lebesgue lemma. 
From Eqs. and pS]). we obtain 

± = j d 4 x(T 12 (x)T 12 (0)) ~±Y, C ^ T ^°»° 

= ^(Ti 2 2 >o-i(ro„)o-^(F 2 )o, (25) 

in the leading order perturbation theory. The right hand 
side of Eq. ([23)) includes only lattice observables. There- 
fore, we can evaluate rj/r^ directly by lattice simulations. 

V. LATTICE MEASUREMENTS 

A. Formulation on the lattice 

In this section, we analyze the ratio rj/r^ on the lattice 
by evaluating three observables on the right hand side of 
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Eq. (|25|) . We have performed the lattice simulations for 
SU(3) pure gauge theory with the standard Wilson gauge 
action. In Table U we list lattice parameters used in this 
work. To investigate the lattice spacing and spatial vol- 
ume dependences, simulations have been carried out on 
four isotropic lattices of volume N%. The lattice spac- 
ing, a, in Table U for each j3 is determined by the string 
tension yfa = 460MeV and the parametrization of ay/a 
in Ref. [20]. All observables discussed in the next sub- 
section, however, are dimensionless and do not depend 
on this physical scale. We use T c / 'yfa = 0.63 [27|, Hi] 
and normalize T by T c . The temporal length aN T = l/T 
of each lattice corresponds to the range of temperature 
0.5T C < T < 4T C , which covers those realized in heavy ion 
collisions at RHIC and LHC. Gauge configurations were 
updated by heatbath and overrelaxation algorithms. We 
have generated 300,000 — 2,000,000 configurations for 
each parameter. Statistical errors were estimated by the 
jackknife method with bin sizes in the range of 50— 1000. 

In the lattice calculations, we employed the traceless 
operator for the energy-momentum tensor, 



T, 



2Tr 



Fp,pF U p , & F pa F pa 



(26) 



For the field strength F^ v on the lattice, we have used 
the clover type plaquette. With this definition, Eq. (|25|) 
reads 



V_ 



-^(^12)0 



l(T v)o-l(F 2 ) 



(27) 



where we have replaced the traceless operator (Too)o with 
— I (e + P) . The definition Eq. (j2"B)) is sufficient to eval- 
uate ti/tk as explained below. To perform the vacuum 
subtraction, we regarded the lattice with the largest N T 
for each /3 to be the vacuum one and subtracted the ex- 
pectation value on this lattice from the one at each N T . 
In fact, the lattice with the largest N T for each f3 corre- 
sponds to a temperature well below T c , where medium 
effects on the expectation value are expected to be well 
suppressed. 

We also used the energy-momentum tensor defined 
above to evaluate e + P. In the calculation, simply their 
expectation values — (Too)o + (133)0 were taken. We note 
that the effect of the trace anomaly is canceled out in this 
observable, and hence the traceless definition in Eq. ([26]) 
is enough for this measurement. While our definition 
for e + P is different from that used in the conventional 
analysis of lattice thermodynamics, the difference should 
be of higher order in lattice spacing and converge in the 
continuum limit. 

Finally, we have measured the gluon condensate (F 2 ) 
from the gauge action using the relation, 



(l?ll= N ± T ^S£ = -6aN; 
T 4 dT 



dg 



da 



(28) 



v 
> 




(b) 



Q_ 

+ 



1 

0.5 



-0.2 
-0.4 



(c) 



p=6.499, L„=1.5fm — e- 
P=6.205, L =2.2fm — ^- 
p=6.000, L_=3.0fm 
p=6.000, L CT =1.5fm — ^- 



0.5 



1.5 



2 2.5 
T/T„ 



3.5 



FIG. 1: (a) Kubo's canonical correlation (T^ 2 ) with only vac- 
uum subtraction. Lattice spacing and volume dependences 
are small, (b) (e + P)/4 obtained from direct measurements 
of (Too)o and (T33)o defined by the field strength constructed 
from clover type plaquette. The results qualitatively repro- 
duce well-known behavior of previous lattice works, (c) Gluon 
condensate measured by conventional approach with standard 
plaquette. 



where Sq is the lattice gauge action constructed from the 
standard plaquette. For the values of adg~ 2 / da, we have 
used those determined in Ref. 123. 



B. Numerical results 

In Fig. [U we show the temperature dependences of 
each term on the right hand side of Eq. (|2"T1) . (T^o , 
(e + P)/4, and (F 2 ) /6. Fig.QJa) shows that the canoni- 
cal correlation {T 2 2 ) a takes negative values for (3 = 6.000 
at T < 2T C . While the negative values are also observed 
for N T — 4 on finer lattices, for N T > 6 the values of 
(T 2 2 )o are consistent with zero within statistics. The er- 
ror bars in Fig. (Ua), however, can give a constraint on 
the upper limit on the absolute values of (T 2 2 )o- Thermo- 
dynamic quantities presented in the lower two panels of 
Fig. Q] show that statistics used in the calculation is suffi- 
cient for their determination. We have checked that nu- 
merical results in Fig. [He) reproduces the trace anomaly 
measured by the conventional procedure well (see, for 
instance, Ref. |27]). The spatial volume dependence of 
these results is estimated by comparing numerical results 
on the two lattices with f3 — 6.000. Fig.[T]shows that the 
spatial volume dependence is almost negligible in our re- 
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2.5 
2 

1.5 
1 

0.5 



P=6.499, L_=1 .5fm 
P=6.205, L c =2.2fm 
P=6.000, L n =3.0fm 
P=6.000, L a =1 .5fm 



0.5 1 



1.5 



2.5 



3.5 



T/T r 



FIG. 2: Temperature dependence of the shear viscosity to 
the relaxation time ratio. The vacuum contribution and the 
contact terms are subtracted at each temperature. A rapid 
increase is seen around the critical temperature. 



0.75 



0.5 



0.25 



=6.499, L n =1 .5fm 
=6.205, L„=2.2fm 



(3=6.000, L n =3.0fm 
3=6.000, L a =1 .5fm 



0.5 1 



1.5 



2.5 



3.5 



T/T. 



FIG. 3: Temperature dependence of the square of the char- 
acteristic speed v 2 — ^/-^(e + P). The results show that 
v 2 ~ 0.25 for the range of temperature T > 1.5T C . 



suits except in the vicinity of T c . 

In Fig. f2J we show the temperature dependence of the 
ratio ry/rjr given by Eq. (|2~T1) . The figure shows that the 
results with all sets of configurations coincide within the 
statistical error, which indicates that the lattice spacing 
dependence of the present results is small. Comparing 
the two results with /3 = 6.000, we observe that the spa- 
tial volume dependence is also small. It seems that rj/r^ 
has a tendency to increase rapidly in the vicinity of the 
critical temperature, but it is hard to conclude because of 
the statistics, rj/r^ is related to the characteristic speed v 
defined by Eq. ([8]) . We show the temperature dependence 
of v 2 in Fig. [31 v 2 is constantly around one-quarter of the 
square of the light speed, v 2 ~ 0.25, for 1.5T C < T < AT C . 
This result indicates that causality is satisfied in the sec- 



ond order hydrodynamics for this mode in this range of 
temperature. The statistical error, however, increases as 
the temperature is lowered to T c . 

It is of interest to compare the present result on v 2 with 
other theoretical predictions. Grad 14 moment approx- 
imation predicts v 2 ~ 1/6 in the high energy limit [l9| . 
and the AdS/CFT correspondence v 2 = 1/2(2 - In 2) ~ 
0.36 [29| . The leading order result of the projection op- 
erator method provides v 2 = P/(e + P) ~ 0.25 in the 
high temperature region [l8j]. The measured values of 
the characteristic speed are close to the value of the pro- 
jection operator method especially above 2.5T C . 

Finally, let us consider the implication of our results 
on the magnitude of the relaxation time r v . It has been 
conjectured that the ratio of the shear viscosity to the 
entropy density rj/s has the universal lower bound 1/47T, 
which is the value obtained by the AdS/CFT correspon- 
dence [30]. On the other hand, it has been suggested by 
the comparison between the results of dissipative hydro- 
dynamics and the experiments at RHIC and LHC that 
the ratio is less than ~ 5x the lower bound above T c 
(H . If one substitutes these values in Eq. §5§ and use 
v 2 = 0.25, the relaxation time r w is roughly estimated to 
be tv ~ O(10~ 2 ) — O(10 _1 )fm in the temperature range. 
This estimate is almost ten times smaller than the result 
for tv in the pion gas below T c [20l | . 



VI. CONCLUSION 

In this paper, we have evaluated the ratio of the shear 
viscosity to the relaxation time of the shear flux rj/r^ by 
measuring Kubo's canonical correlation of a spatial com- 
ponent of the energy-momentum tensor for the range of 
temperature 0.5 < T/T c < 4 with SU(3) lattice gauge 
simulation. The canonical correlation diverges in the ul- 
traviolet region and we regularized it by subtracting the 
vacuum contribution. There, however, exists an extra 
contribution from the contact terms that remain even 
after the vacuum subtraction due to their temperature 
dependence. Using the leading order OPE result for the 
contact terms and measuring the expectation values of 
the operators on the lattice, we determined their tem- 
perature dependence. Then, we have obtained the tem- 
perature dependence of rj/r^. The number of the free 
parameters in relativistic dissipative hydrodynamic mod- 
els was reduced by this analysis. We also investigated 
the temperature dependence of the characteristic speed 
of the transverse plane wave in the same temperature 
range. We have found that the characteristic speed in 
gluon media is v ~ 0.5 for T > 1.5T C . The fact that the 
characteristic speed does not exceed the speed of light 
within statistical error indicates that relativistic hydro- 
dynamic models based on IS theory retain the causality 
for the transverse modes in the quark-gluon plasma pro- 
duced in heavy ion collisions at RHIC and LHC. Note, 
however, that our calculations are in the quenched ap- 
proximation, i.e. with no dynamical quarks. 
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The approach in this work is also applicable to the 
bulk channel. By repeating the same procedure as in 
Sections 3 and 4, one can obtain the ratio of the bulk 
viscosity to the relaxation time of the bulk flux, C/ T n, 
from Kubo's canonical correlation of diagonal compo- 
nents of the energy-momentum tensor. This is currently 
in progress. 
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